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O ' Abstract 
(N 

The 3D Ising-like system in the external field is described using the non- 
perturbative collective variables method. The universal as well as nonuniver- 
. sal system characteristics are obtained within the framework of this approach. 

! The calculations are carried out on the microscopic level starting from the 

Hamiltonian. They are valid in the whole h — T plane of the critical region. It 
q ■ is established, that the contributions related with wave vector values k — ► 

exhibit the properties of the total system near the critical point. The be- 
haviour of the susceptibility as function of the temperature in the presence of 
the field is investigated. The locations of the maximums susceptibility on the 
temperature scale for different values of the field are established. 
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& ' 1 Introduction 



Despite the simplicity and clearness of the 3D Ising model, it is still the ob- 
ject investigated with help of the approximate methods jTj. We propose the 
description of the behaviour of the Ising-like magnet in the external field by 
collective variables (CV) method (2j. This approach is similar to the well- 
known non-perturbative Wilson-Polchinskii Renormalization Group method 
[HIE]. In particular, the coarse-grained procedure with using renormalization 
group (RG) transformation is employed in the both methods. But, in contrast 
to the Wilson approach pj, the CV method does not involve any phenomeno- 
logical assumptions and adjusting parameters. It is based on the physically 
well-grounded and mathematically rigorous use of the space of CV, which are 
associated with modes of spin moment density oscillations. As result, the 
method allows one to obtain the macroscopic characteristics as functions of 
the microscopic parameters, contained in the system Hamiltonian. 

In the 3D spin systems exhibiting the critical behaviour, the presence of the 
asymmetry (presence of an external field) sufficiently complicates the problem. 
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Usually, the series in scaling variable, which are some ratio of the field and 
temperature variable, are used for description of the system behaviour. In 
this case, one of the factors dominates and forms the system critical behaviour. 
Using CV method, we have already obtained such characteristics in the form of 
a power series in the scaling variable in the regions of the weak and strong fields 
for temperatures above and below T c (T c is the phase transition temperature 
in the absence of an external field) (HI Ej. When the role of the field and 
temperature is equally important, the series are not valid. For complete system 
description, it is necessary to perform calculations in the whole h — T plane, 
where T is the temperature, and h is an external field. 

In the perturbative approaches (e-expansion, high-temperature expansion, 
etc), the parametric equation of state is used for description such systems [HJ- 
Within the framework of the Wilson-Polchynskii method, the investigation of 
the 3D Ising model critical behavior at the presence of an external field is per- 
formed in the so-called local potential approximation (LP A) In particular, 
the critical exponents of the correlation length v and order parameter S as 
well as asymmetric leading corrections to scaling related with presence of the 
field are calculated. The description of such a system can be used also for in- 
vestigations of the nonmagnetic systems. For instance, the critical behaviours 
of the system "gas-liquid" in the plane (temperature,pressure) pQ, Standard 
Model of the electroweak interactions in the plane (higgs mass, temperature) 
^Ol^B) were shown to have similar features of the phase transition. These 
models belong to the 3D Ising universality class. 

2 Basic relations 

We represent the calculations, which are valid in the whole h — T plane of the 
critical region. Starting with the Hamiltonian 

H = ~J2^( r ^i- h J2^ (i) 
z y i 

and using the "layer-by-layer" integration of the partition function, we obtain 
the critical exponents, free energy and other thermodynamic and structural 
characteristics of the system. The quantities a\ are z-components of the spin 
operator, whose eigenvalues take on the values ±1, and h is an external field. 

The interaction between spins on the simple cubic lattice is described by 
the exponentially decreasing function $(nj) = Aexp(— ry/6) of the distance 
Hj between particles at sites j and 1. Such an approach makes it possible to 
investigate the system behaviour depending on the microscopic parameters, in 
particular, on the radius of the effective interaction b and constant A. 

For the following calculations, it is necessary to perform the transition from 
quantities a\ to the CV pk through the functional representation for operators 
of the spin density oscillation modes 

j o k = (ViV)- 1 ^ f 7 1 e X pHkl). 
l 

Here N is the number of particles in the system. The detailed procedure 
of this transition is described in works |2j. We use also the parabolic 
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approximation of the Fourier transform of the interaction potential in the 
form 

1J ~ \ $ = *(0)$, B < k < B, (Z) 

where B = tt/c is the boundary of the Brillouin half-zone, and the quantity 
Bo = B/sq and so > 2 is the starting division parameter. This parameter 
defines the boundary of the region, where the parabolic approximation of the 
interaction potential is valid. The quantity $ is the constant and <&(0) = 
8-7rA(b/c) 3 , where c is the lattice constant ^3j- We begin the calculation with 
partition function EH 

Z oc / (dp)% (dcv)% e W h J2 (*(*) " $oW-k 
J L Z kGBo 

x exp(2vri ^ p k u> k )Z(u>). (3) 
Here /3 = and k is the Boltzmann constant. The quantity 

z(u) = ^Jy[ { ~ 2 ^ )n N i - n ' 2 



-M n (ti) ^k 1 ---Wk,A 1 +...+kJ \ 



(4) 



is the the cumulant expansion of the Jacobian of transition from the set of 
the spin variables a\ to the set of CV /5 k , where <5 kl +...+ kn is the Kronecker 
symbol. The coefficients M. n are defined with help of the relations: 

an \ a 

A^ = lncosh/ l / , M n = -^. (5) 

The integer even values of the quantity n in formulas (@J and (jSJ determine 
the choice of the effective potential (models p 4 , p 6 , etc), which will be used for 
the description of the system. The value n = corresponds to the mean-field 
approximation, and at n = 2, we have Gauss distribution of the fluctuations 
in the system. The odd degrees of CV appear as result of the presence of an 
external field. In this work, the calculations are carried out employing the 
simplest non-Gaussian model with n = 4. The results allow one to take into 
account long-range interactions, which determine the critical behaviour of the 
system. It should be noted, that the variable po, which corresponds to the zero 
value of the wave vector, is crucial in the forming of the partition function 
© . The average value of this quantity is the order parameter of the system 

m 

For the integration of ((HI, one should pass to the following field variables 
by performing the substitution: 



Pi = (V^Vo) 1 p k exp(ikl), 
keBo 

= (\//Vo) _1 ^kexp(-ikl), 
keBo 
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where iVr 



o 



NSn 



In this case, the partition function takes on the form 



Zoh J (dp)o (cIuj)q exp 



~(3 (*(*) ~ $o W-k 



keB 



(6) 



The vector 1 varies in the volume of periodicity V = NqCq, where cq = sqc, 
and 







2 W exp [iV(-^ + Mo)], 
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It should be noted, that the variable uj\ can be interpreted as the scalar field, 
which is similar to the Habbard-Stratanovich field. Within the frames of the 
quartic model, the integrals I\(pi) is defined by the expression 



h(f>i 



l 



irV2 



du!\ exp[27ripi( 



d/2,, 

2vr 



)] 



r 2 " 4i 
x exp[wi - C h UJ y \. 



(7) 



Here a h = s^h /2 and c h = Sq (1 — 4/i ), We perform also the substitution of 
variables to cancel the cubic term in the cumulant expansion. The coefficients 
in the exponent expression are approximated by power series in field variable 
to h' 2 . This approximation is quit precise, since in the vicinity of the critical 
point, the maximum values of the field are of order of 0.01. Such values of 
the field correspond to the real fields, which are commensurate with fields of 
the order parameter saturation of the ferromagnetic materials. We will find 
the result of the integration of the expression Q in the form 



e a ° exp 



n=l 



(8) 



where the coefficients a n are defined using the condition 



d n K x 



Pi=0 



dpf 



Pi=0 



Since the quantity cj| < 1, we present the integrand of the Q by the power 
series in the quartic term. Such an approximation is used only in the space of 
variables p\. For the coefficients a n , we have 

a = ln(>/7r(l - 



a i 



-4 /2 h', 



a 2 = l- 3c h , 
a 3 = 0, a 4 = 6c h . (9) 

In the space of the variables pk, the partition function takes on the form: 

Z = Z j (7TV2)- 1 e^ +a ^ No J(dp) N ° expf-axViVopo 
1 



- d(k)p k p_ k 



kGBo 

"4!iVi 



+k 4 \ 



(10) 



keB 
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where d(k) = a 2 + /?<3?o ~~ /3$(k). The presence of the external field is repre- 
sented only by the linear term in the exponent of the expression (jlT)]) , It should 
be stressed that the same result is valid for the models with n = 6, 8, .... The 
field dependence of other coefficients is very weak in comparison to coefficient 
in the linear term. 

The partition function lfTT)|) is starting point to perform the "layer-by-layer" 
integration with using RG transformations for the quantities a n . Beginning 
the integration from the variables p k with large values of k, we divide the the 
phase space of the CV p k into the layers with the division parameter s. In 
each nth layer (where B n+ \ < k < B n , B n+ \ = B n /s) the Fourier transform 
of the potential <3?(A;) is replaced by its average value. Such an approximation 
corresponds to the LP A in the Wilson-Polchinskii nonperturbative theory [3j. 
The detail description of this procedure is given in the paper ^Hl- The result 
of the integration can be represented by the expression 

Z = Z Q Q 1 ...Q n [Q(P^)] N ^I n+1 , (11) 
where Q = [(■KV2)- l e- a 'h +a °Q(d)] N ° , and quantities 

Qn = [Q(P {n - iy )Q(d n )] Nn (12) 
are the partial partition functions of the nth block lattice. Here 
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Q(p(»-i)) = (2^P 2 (n - 1) )-5(l - -G^), 



Q(4) = (^y)"7l(l-74 n) ) (13) 



as well as 
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H"~ 1} = -At 7(1 + ^4"-" 



Sn-l) 

G (n-i) =s -dG (l + G 2 h% l ~ 1) ), 

(4 W) )V 2 

The quantities 7, 71, Go, t 2 and G 2 do not depend on the field. They are 
given in The coefficient 

d n (B n ,B n+1 ) =d n (0)+qs- 2n 

is defined through the value, which related to the interaction potential f3Q(B n , B n ^ 
averaged on the interval k £ B n \B n+ i. As result, we have 

d n {k) = 4 n) + 0$(B n , B n+1 ) - 0$(k). (15) 

Here q = qP^(0) and q = (bn / c) 2 Sq 2 (1 + s~ 2 ). The integral I n +i in the 
formula ifl!]! has the following form: 



In+i = f(dp)% +1 expl -of +1) \/iV n+ ip 



(n+l) 

H 
4!iV, 



a (n+1) 1 

iw— £ p kl ...p- k4 4 1+ ... +k4 |. 

KGCn+l / 



(16) 



1; 
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The coefficients af"^ are defined by the recurrent relations (RR) 

4 n+1) = /oo(4 n) ) 1 / 2 (i + a 2 4 n) ) 

al n+1) = S - d / iaf(l + a 4 4 n) ), (17) 

obtained as result of the RG transformations. The quantities a 2 , «4, /oo and 
/oi are given in [13J. Due to the presence of the field, we obtain I + 1 RR for 
the models p 21 , where I = 1,2,3... Performing the transformation p^ = sp'^, 
introducing the notation oj n = s n a[ \ r n = s 2n d n (0), u n = s An a^ and taking 
into account the expression for h 2 1 ^ in (dJ, the equation lfT7|l can be written 
in the matrix form 

' Rn \ / uj n - u* 
R22 R24 r n -r* 
R42 Ru J \ u n - u* 

where R n = s^r , R 22 = s 2 f 00 a 2 V6, R24 = s 2 f 00 u*- 1 ^, R 42 = sfo^VEu* 1 / 2 
and -R44 = s/oi- The quantities w* = 0, r* = —q and u* = q(l — s~ 2 )//oo are 
the fixed point coordinates. For parameter s, there exists a preferred value 
s = s*, which is determined by condition = = 0. 




3 RR Solutions. Scaling Region. 

The fixed point is associated with the phase transition point, since in this case, 
the fluctuations become infinitely large, and it is necessary to perform infinite 
number of iterations for system description. But at the fixed values of the 
field h! and reduced temperature r = (T — T c )/T c , this number become finite. 
Thus, one can linearize the RR (|T7|) in the case of the small deviations from 
the fixed point and find the coefficients of the nth block structure through 
their initial values Q- Hence, one can find the eigenvalues and eigenvectors 
of the RG matrix and write down the solutions of the equation lfT%|) in the 
following form: 



r n : 

Un 



UJ 



J/2,, 



h'E^ 



1/2/ 



-1 jpn 
4 > 



2^(0), 1/2, 



(19) 



Here E\ = 20.98, E2 = 7.37, E4 = 0.4 are the eigenvalues of the RG matrix, 



l ' 2 R 2A {E A - R22)- 1 , T$ = u*- 1 / 2 R 42 (E 2 - Ru)' 1 - The coefficients 



-"24 

cj^ and 0^2 are expressed by relations 



JO) 

-kl 



Ck2 



l-q>-q- T 2 { fu <p 1/Z (/?<D(0)/3 C <1>(0)) 



1/2/ 



_ T (0) 1/2' 
J 24 



r (0) T (0) 
-24 1 42 



no-n*-T 4 ( 2 °Vo /2 ^(0)(ro-r*; 



T (0) T (0) 
1 24 1 42 



6 



where tpo = (q(l — s* _2 )//oo) 2 - These values do not depend on the field. 

Using the linearity condition for relations (jlH]) , we define the number n of 
the iterations. This is so-called scaling region, which is defined by the values 
of the field h! and temperature r. The deviations from the fixed point are 
formed mainly in the first two equations of the system ifT^l. since in the last 
equation for u n , the coefficient near the quantity tE 2 is small in comparison 
to the coefficient near the same quantity in the equation for r n . Taking into 
account above mentioned condition, we write down the equation 

(-s d /2 h'E r ^ +1 ) 2 + (cgV/^oK^ 1 ) 2 = r* 2 , (20) 

which allows us to find the number n p of the "layer-by-layer" integrations. 
The definition of the quantity n p = n p (r, hi) as function of the field and 
temperature is essential for investigating the crossover between temperature- 
dependent and field-dependent critical behaviour of the system. Indeed, this 
number determines the size of the block structure, which is commensurate with 
correlation length. Equality of the terms in the left side of the equation (|2T)|) 
corresponds to the state of the system in the vicinity of the line determined 
by relation 

roc/i^, (21) 

where (3 and 5 are the critical exponents of the order parameter. In this 
case, the field and temperature equally influence on the forming of the critical 
behaviour of the system. 



4 Universal and Nonuniversal characteris- 
tics of System 

In the high-temperature region, the partition function is represented as 

Z = Z QoQl...Qn p Qn p+ l[Q( J P (n+1) )]^ +1 /n p+ 2. (22) 

It should be noted, that in comparison to (0, the expression (|2"2|l contains 
additional partial partition function Q n +\. Its presence is related with per- 
forming the additional step of integration. Such an integration is carried out 
for reaching the dominant value of the coefficient d np+ 2 in comparison to other 
coefficients in I np +2 for any k ^ 0, where k G B n+2 (see iflfijl ). Thus, for these 
values of the wave vector, we perform the integration with respect to the 
collective variables in the Gauss approximation. 

Using the partition function l[22j), we represent the system free energy as 

F = F^ + F k=0 . (23) 

The second term in this expression is related to the integration of the quantity 

= /(^o)exp N l / 2 h' PQ - l -d rip+2 {Q)pl 



(24) 



with respect to the variable po using the steepest-descent method. As was 
mentioned above, this variable corresponds to the order parameter and forms 



7 



the main contribution to the free energy of the system. Using the substitution 
of variables po = y/N po and the maximum condition for the expression in the 
exponential function (see 1|24|)) by differentiating with respect to the variable 
p~o, we find the root of the cubic equation in the form 

np + l 

Po = &os~ 2 • (25) 

The choice of the valid root is based on the free energy minimum condition. 
The term i*k=o is defined by expression 



(26) 



The expression l|25|) for p$ describes the order parameter behaviour and can 
be interpreted as equation of state. 

For another term in the right side of l|2~3j) . we have the relation 

F k ^ = -NkT [in cosh ti + l + hh 2 

+h f + / 4 f2 + Z W eS -3(n p +l)] j (27) 

which is the result of the integration with respect to the variables related with 
nonzero values of the wave vector. Here, the coefficients Iq , lyr e , 1%, h and 
h are independent of the field ^H], h = h'/q, f = cf^r/q. The contribution 
l|2*7|l is small in comparison to -Fk=o- Indeed, its magnitude depends on the 
accuracy of the integration. With increasing the accuracy of the calculations, 
this contribution vanishes. It is physically grounded, since as result of the 
RG transformations, we deals with larger lattice structure, which exhibits the 
properties of the previous smaller ones. Thus, the result of the last integration 
with respect to the variable po should represent the properties of the whole 
system. In the figure ^ the average value of the variable po and order param- 
eter obtained by differentiating the free energy F with respect to the field as 
functions of the field are represented. As one can see, the difference between 
these quantities is neglible small. The result for the order parameter agrees 
with the data obtained by the Monte-Carlo simulations [16J. Some deviation 
in the strong field region is explained by the insufficient accuracy of the model 
p 4 . As has been shown in the similar investigations in the absence of the field, 
the better results can be obtained with using the higher approximations for 
the effective potential, in particular, p 6 model [T31 fT7| fT%]. 

Thus, the expression (ffiSjl represents the total free energy of the system 
obtained using the sequential mathematical transformations starting from the 
Hamiltonian (QJ. The relation l(23]l is the equation of state, which is valid in 
the whole h-T plane. 

In the case of h = 0, taking into account the solution of ifflljl for n p , the 
free energy takes on the form: 

F k=0 = —NkTRf 3 ", (28) 

where R is the sum of the last two terms in (J2(ij). The critical exponent of the 
correlation length satisfies the relation 

In s* , . 

— E^-Ml. (29) 
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The equation of state is transformed into expression 

/ o = (To f2. (30) 

At T = T c , the second term in the ifflty disappears. In this case, the free 
energy can be written down as 

F k=0 = -NkTR'h'l , (31) 

where R' = a r*/s^ /2 - r* S - 2 a$/2 - sgs" V<rg/4! and h' = h's d Q /2 /r*. Corre- 
spondingly to f|25[l . the equation of state has the form 

p o = a o 0. (32) 

The critical exponent of the order parameter takes on the value 5 = 5. Within 
the frames of the LPA of Wilson-Polchinskii theory, the values of critical expo- 
nents of the correlation length and order parameter are v = 0.6895 and 5 = 5 
respectively [9 J (one of the best estimates for the exponent of the correlation 
length is v = 0.63002 [19J). Using the expression l(26|) for the free energy, 
we calculate also the system susceptibility %■ I n the figure it is shown 
the dependence of this characteristic on the temperature for different values 
of the field. The maximums of sussepibility are located at the value of the 
scaling variable a m = T m /h l ^ & as 0.77. The same quantity in the mean-field 
approximation takes on the value as 0.8255 |2()j . 



5 Conclusions 

The general description of the 3D Ising-like system in the external field by 
the CV method is represented. The approach is nonperturbative and can be 
used for the various systems, for which the perturbative methods are not quit 
appropriated. It allows one to calculate macroscopic characteristics of the sys- 
tem derived from the microscopic quantities. The investigations are performed 
in the whole h-T plane of the critical region. The main advantage of the CV 
approach in comparison to others is the absence of any phenomenological as- 
sumptions and adjusting parameters. All calculations are carried out on the 
microscopic level. 

Using the coarse-graining procedure, we obtain the explicit analytic expres- 
sions for the free energy and order parameter of the system. It is established, 
that the contribution related to the po variable exhibits the properties of the 
total system. Another part corresponding to the nonzero values of the wave 
vector is essentially small and depends on the accuracy of the integration. 
The crossover to the limit cases {h = 0, T ^ T c ) and {h ^ 0, T = T c ) is 
demonstrated. The critical exponents of the order parameter and correlation 
length are calculated. The dependence of the order parameter as function 
of the external field is plotted. The result agrees with the result, obtained 
by Monte-Carlo simulations, in the weak-field and crossover region. Some 
deviations in the region of the strong fields can be explained by the quartic 
approximation of the effective potential. In order to obtain more accurate 
results, it is necessary to employ the approximations with n = 6,8,... The 
system susceptibility is calculated and plotted as function of the temperature 
for different values of the field. We investigate the maximum location of this 
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function. It is shifted to the region of the stronger fields in comparison to 
result obtained using the mean-field approximation |2()| . 

The investigations can be extended to the model with the n-component 
order parameter as it was shown in in the case of the absence of the 
field. The proposed approach permits one to study also the systems with limit 
number of the particles. 
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0-6 - — ■— Monle-Carlo result [13] 
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Figure 1: The order parameter of the system as function of the external field at 
r = 5-10" 3 . 
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Figure 2: The system susceptibility as function of the temperature for different 
values of the field. 
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